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Abstract. A complete phase diagram of Zn substituted quantum quasi- two- 
dimensional helimagnet LiCu202 has been presented. Helical ordering transition 
temperature (Th) of the original LiCu202 follows finite size scaling for less than 
^ 5.5% Zn substitution, which implies the existence of finite helimagnetic domains 
with domain boundaries formed with nearly isolated spins. Higher Zn substitution > 
5.5% quenches the long-range helical ordering and introduces an intriguing Zn level 
dependent magnetic phase transition with slight thermal hysteresis and a universal 
quadratic field dependence for Tc{z < 0.055,H). The magnetic coupling constants of 
nearest- neighbor (nn) Ji and next-nearest-neighbor (nnn) J2 {(y=J2/Ji) are extracted 
from high temperature series expansion (HTSE) fitting and N=16 finite chain exact 
diagonalization simulation. We have also provided evidence of direct correlation 
between long-range helical spin ordering and the magnitude of electric polarization 
in this spin driven multiferroic material. 



PACS numbers: 75.10.Pq, 75.40.-s, 75.50.Lk, 75.60.-d, 75.85.+t, 76.60.-k, 77.22.Ej 
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1. Introduction 

LiCu202 is a complex spin-driven multiferroics which shows spontaneous electric 
polarization below the spin spiral ordering temperature near 22 K.[Lj The noncollinear 
spiral spin structure allows inversion symmetry breaking and manifested a finite electric 
polarization. The nature of helimagnetic ordering and its correlation to the existence 
of electric polarization has been the central topic in the research of multiferroics. In 
particular, microscopic model based on spin supercurrent or inverse Dzyaloshinskii- 
Moriya interaction provided the theoretical ground for the generation of electric 
polarization P within noncollinear magnets, i.e., P oc Q x (S^ x Sj), where Q is 
the cycloidal vector pointed from neighboring spin sites i to j.O [3] Incommensurate 
helimagnetic ordering has been found in LiCu202 and accompanies weak electric 
polarization of multiferroic nature, although confiicting results on the cycloidal spin 
plane assignment and disagreement between experimental ferroelectric properties and 
spin-current model predictions remain. [H [3l HJ |6l [7j 

LiCu202 has a Cu-0 chain structure formed with edge-sharing CUO4 plaquettes 
in the ab plane while these chains are connected through CUO2 dumbbells along the 
c direction as shown in figure [l] [8j This compound is uniquely composed of nearly 
equal amount of Cu+ and Cu^+, with Cu+ sitting in the 0-Cu-O dumbbell (along the c 
direction) to connect the nearest-neighbor edge-sharing spin chains. Undoped LiCu202 
shows a helical spin ordering below ^ 22 K, as indicated by the dx/dT peaks, which 
has been verified by magnetic neutron scattering studies before. [HI [4] Two dx/dT peaks 
occur near ^ 22 and 24 K, which are anisotropic and more pronounced for magnetic 
field applied || ab and || c, respectively. The anisotropic magnetic correlation lengths ^ab 
and have been extracted using synchrotron resonant soft x-ray magnetic scattering. 
It is suggested that the transitions near ^ 24 and 22 K correspond to the quasi- 
two-dimensional short-range and the following three-dimensional long-range orderings 
respectively. [21 [9] Li deficiency is proposed to introduce spin- 1/2 Cu^+ to the bridging 
spinless Cu+ site and the helimagnetic ordering temperature is reduced only slightly 
from ^ 22 K to 20 K for samples nearly free of Li deficiency. [TO] 

Many early works attempted to describe LiCu202 as an alternating Heisenberg 
antiferromagnetic chain and applied Bonner-Fisher fitting to the magnetic susceptibility 
for T > Th] however, the fitting is generally unsatisfactory. [TTl [T2] Later nearest- 
neighbor exchange coupling Ji is found to be ferromagnetic and J2 to be 
antiferromagnetic, based both on LDA calculation and neutron scattering spin wave 
analysis, which satisfies the classical spin spiral ordering requirement of IJ2/J1I > | as 
verified experimentally. [131 El Quantum fiuctuation must sustain in the spin spiral 
ordered state to account for the discrepancy found in polarized neutron scattering 
intensity for this frustrated low dimensional spin 1/2 system. [3] 

If we assign the shortest Cu-Cu distance within the edge-sharing chain to be the Ji, 
its Cu-O-Cu bond angle of ^ 94° is expected to be ferromagnetic interaction, following 
the Goodenough-Kanamori- Anderson rule, as summarized in figure [2| [TSi [TBI ITT] There 




Figure 1. (Color online) (a) Crystal structure of LiCu202. (b) Magnetic couplings 
described by linear chain model and zigzag chain model. The purple dot line represents 
the upper chain within layer A. The orange dashed line represents the lower chain 
within bilayer B. Linear chain model is adapted from Gippius et al. [13j and zigzag 
chain model is adapted from Masuda et al. [14 distinguished with superscript M. Note 
the Jw is defined as superexchange coupling between upper and lower Cu^+ in bilayer 
A bridged by the zigzag shape channel through Cu+. 



have been two major views on assigning spin chains and its interchain coupling before, 
one simply viewing the edge-shared CUO4 network as linear spin chains shown in 
figure [l](b),[Il [IB] and the other connecting the Cu between linear CUO2 chains to 
form zigzag chain described by the double layer-B as shown in figure [l](b) . jlj Zigzag 
chain description is only possible when non-negligible coupling J^i defined in figure 
[l](b) exists; [4J however, LDA calculations and alternative neutron scattering spin wave 
analysis suggest it to be negligibly small. [131 E] the other hand, based on the 
model in figure [l](b) , the interchain coupling between linear chains through the Cu- 
0-Cu superexchange has been estimated to be comparable to J2 based on LDA 
calculations. In particular, Wannier function modified Hubbard model suggests strong 
orbital overlap between oxygen and the bridged Cu+ ions, as the Jw defined in figure 
[l](b).[19] These two kinds of structural views led to different assignment of Ji, J2 
and interchain coupling J^, but the discrepancy has been resolved based on neutron 
scattering spin wave analysis of alternative fittings. [T^ Herein we use the linear chain 
model with non-negligible interchain coupling in the following description and define 
Ji, J2, c^=J2/Ji and interchain coupling J± based on the model shown in figure [l](b). 
When J± is non-negligible as predicted from both calculation and experiment, J± (or Jw) 
would connect linear chains equally from both sides along the a direction and convert the 
system from one-dimensional to quasi- two-dimensional. In fact, such two-dimensional 
view has been confirmed from its renormalized classical behavior of short-range in- 
plane magnetic correlation length based on resonant soft x-ray scattering measurement 
results. j9] 
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Figure 2. (Color online) Nearest neighbor coupling Ji versus Cu-O-Cu bond angles 
summarized for some samples with edge-sharing Cu02 chains, adapted from references 
[inilTT]. The inset shows values up to 180° for mostly two-dimensional layered structure 
and the main plot amplifies data near 90° region for the chain compounds. All Ji's are 
obtained from neutron scattering experimental studies, except those from susceptibility 
fitting and labeled as empty square in blue. 



Isotropic spin- 1/2 Heisenberg antiferromagnetic chain has been the most widely 
studied model since the introduction of the Bethe ansatz in 1931. [20] There have been 
considerable theoretical works on the spin- 1/2 chain system with both nearest-neighbor 
(nn) and next-nearest-neighbor (nnn) of antiferromagnetic interactions as described by 
the Hamiltonian of 

H = Ji ^(s^ • s^+i) + J2 ^{^i ' Si+2). (1) 

However, studies on chains with ferromagnetic Ji and antiferromagnetic J2 interactions 
have been rare, and it was only until recently that great interest was re-focused on 
the quantum spin system with emerging electric polarization concomitant of helical 
spin ordering, such as LiCu202, LiCuV04 and Li2ZrCu04 as summarized in figure 
[2|yj, m [211 [22] Helimagnetic spin ordering for a spin chain occurs when frustrating 
Ji and J2 satisfy a unique condition of IJ2/J1I > | classically. We may also 
view spin periodicity increases from two (for antiferromagnetic ordering) to infinity 
(for ferromagnetic ordering) with helical ordering of intermediate commensurate or 
incommensurate periodicity in between in a classical picture, as described in figure 
[3] by its pitch angle 9 = arccos(- Ji/^ J2) growing from zero ( J2/-^i=-l/4) to tt {J2/Ji= 
1/4). On the other hand, finite chain calculation on one-dimensional quantum spin-1/2 
Heisenberg system suggests the occurrence of a gapped two fold spin singlet (dimer) 
ground state at J^^^/J^^=l/2 (called Majumdar-Ghosh, MG, point) and a first order 
phase transition at J2/J1 = -1/4 (called FF point) before entering the ferromagnetic 
regime. [23j The critical ratio of J^^ /J^^ to have a gapped phase from the gapless 
"spin liquid" (^2=0) state for the spin 1/2 antiferromagnetic Heisenberg chain system 
has been estimated numerically to be ^0.2411. [24] Finite chain numerical calculations 
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Figure 3. (Color online) T = phase diagram of one-dimensional quantum spin 1/2 
Heisenberg system with nearest-neighbor Ji and next-nearest-neighbor J2 isotropic 
exchange constants. The Hamiltonian is defined in equation [l] Figure is reproduced 
following Bursill et al\2Z\ 

also suggest the occurrence of complex periodicity when approaching the FF point. [23j 
There are very few quantum spin chain systems that sit within the second quadrant 
of the phase diagram shown in figure |3| except that the Li2ZrCu04 is claimed to be 
the one closest to the ferromagnetic critical point. f22l [25j| LiCu202 of incommensurate 
helical spin ordering with confirmed ferromagnetic Ji could be the best candidate so far 
to explore this region. [21j 

Nonmagnetic ion substitution into the low dimensional quantum spin system 
has been proved to be not just a simple dilution eflFect. For example, analytic and 
quantum Monte Carlo studies on nonmagnetic ion substitution to the Cu site in 
two-dimensional antiferromagnet La2Cu04 is shown to induce substantial frustrating 
interactions, approximately > 0.4 J per impurity, and can resolve the discrepancy 
between the experimental and fittings by previous theoretical models. [2^ More 
interestingly, Zn substitution of Cu in the quantum spin chain CuCeOa induces an 
unusual antiferromagnetic ordering coexisting with the spin-Peierls (SP) state of spin- 
lattice dimerization.pTl l28] The surprising finding of the coexisting antiferromagnetic 
and SP states (dimerized antiferromagnetic ordering) is diflFerent from the conventional 
uniform antiferromagnetic state. It is shown that staggered moments can be introduced 
and with large spatial inhomogeneity in the ordered moment size (correlation length 
^ 10a) through spin chain perturbation. [291 [30] Clearly nonmagnetic spin zero 
perturbation is an eflFective way to probe the mysterious low dimensional quantum spin 
system. 

Herein, we report another surprising finding on Zn substituted quasi-two- 
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dimensional quantum spin-1/2 system LiCu202. Below ^ 5.5 % Zn substitution, helical 
ordering temperature Th(x) reduces as a function of Zn level and follows an interesting 
finite size power law, which implies an intriguing domain formation. Novel phase 
transition below ^ 20 K is found by > 5.5 % Zn per Cu02 chain substitution. Magnetic 
susceptibility measurement results suggest that a new phase transition emerges with 
a character of spin gap opening with a magnetic field dependence which is similar to 
that found in a spin-Peierls phase transition. We tentatively propose one spin dimer 
model which is consistent to the available magnetic susceptibility, specific heat and 
NMR measurement results, a satisfactory model fitting and simulation is given as well. 
The spin gap size is found to be roughly proportional to the Zn substitution level, and 
these isolated dimers are suggested to sit on a background of frustrated spins in the 
finite spin chains cut short by the Zn defects. 

2. Experimental details 

A complete series of single crystal LiCu2-2Zn202 with z ^ 0-0.10 have been grown 
using traveling solvent fioating zone (FZ) method. The grown crystals with chemical 
analysis are summarized in table [T] Feed rod of nominal Zn levels of 0, 1, 3, 4, 4.5, 
5, 10 and 15% are prepared through solid state reaction route starting from Li2C03 
: CuO : ZnO mixture of molar ratio 1.2 : 4-z : z, where 20 % excess Li is added to 
the initial to compensate for the Li vapor loss. The feed rod is annealed at 750 °C 
after the thoroughly ground powder mixture is treated at 850 °C under O2 fiow for 12 
hours each for one time with intermediate grinding. To reduce further Li loss during 
the molten stage and to prevent Cu^^-rich Li2Cu02 impurity phase formation, high 
pressure 0.64 MPa) argon atmosphere is used during crystal pulling and the pulling 
rate is maintained at 3 mm/hr with 20 rpm feed/seed rods in counter rotation. Li 
content of the Zn-free crystal has been determined using combined thermal gravimetric 
analysis (TGA) and iodometric titration methods as reported earlier. ^lOj and the Cu 
and Zn contents have been determined using combined Inductive Coupled Plasma (ICP) 
and Electron Probe Microanalysis (EPMA). While it is impossible to determine the Li 
content accurately with combined TGA and titration methods for the Zn substituted 
samples, and ICP technique alone cannot provide accurate Li content within 10 % 
error, we assume that Li content does not change under the same high pressure Argon 
atmosphere growth condition and should be maintained near 0.87±0.03 according to 
our previous conclusions on the pure Lia;Cu202 crystal. jTO] Details of FZ growth and 
chemical analysis of Lia^Cu202 single crystal have been reported previously; however, we 
describe Li content to be 1 per formula for simplicity in this paper. The homogeneity 
issue of Zn substitution is examined with EPMA mapping along the growth direction 
within 2 % error within 8 mm and the magnetic phase transition temperature is single 
and with sharp transition width represented by dx(T)/dT peak near phase transition. 

Magnetic susceptibilities are measured with SQUID magnetometer (Quantum 
Design MPMS-XL) with a magnetic field of 1 kOe applied along and perpendicular to 
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Table 1. Zn concentrations determined by EPMA and ICP for LiCu2-z'Z^^z02 crystals 

studied. 

sample ID z (nominal) z (EPMA) z (ICP) 



AnU 


U 


AT A 

IN A 


AT A 

JNA 


Znl 


0.01 


0.016(3) 


0.008(4) 


Zn2 


0.03 


0.034(2) 


0.023(7) 


Zn3 


0.04 


0.038(2) 




Zn4 


0.045 


0.041(1) 




Zn5 


0.045 


0.052(1) 




Zn6 


0.05 


0.058(2) 


0.050(6) 


Zn7 


0.10 


0.081(3) 


0.080(5) 


Zn8 


0.15 


0.109(3) 


0.108(4) 



the ab plane. The helimagnetic spin ordering cannot be identified clearly from magnetic 
susceptibilities, where xC^) drops smoothly at the transition temperature but can only 
be identified clearly by its derivative dx/dT peaks. Part of the x-ray structure data were 
taken using synchrotron x-ray facility NSRRC in Taiwan. For dielectric measurement, 
the crystals were shaped into a thin disk with thickness approximately 50-100 /xm. 
The top and bottom surfaces were then covered by silver epoxy as electrodes. The 
dielectric constant was measured by a high-precision capacitance bridge at 1 kHz. 
At room temperature, the dielectric loss is less than 1,000 nano-Siemens (nS) and 
decreases at lower temperature reaching 0.02 nS below 30 K. The electrical polarization 
is determined by integration of the pyroelectric current which is detected at a rate of 
0.167 K/s after cooling the samples in a poling field of 1200 kV/m from 50 K. ''Li NMR 
measurements were carried out based on standard pulsed NMR techniques. Field-swept 
NMR lineshapes were measured by integrating the spin echo signal at a fixed frequency 
of / = 33.095 MHz while sweeping the magnetic field. Spin-lattice relaxation rate 1/Ti 
was measured by applying an inversion pulse prior to the spin echo sequence with varied 
delay times. 

3. Results and Discussions 

Nonmagnetic Zn ion substitution to the spin 1/2 Cu^+ site within edge-sharing Cu02 
linear chain is expected to disrupt Ji effectively and to cut down the infinite chain 
into even and odd finite chains, although quantum fiuctuation, next nearest neighbor 
J2, and interchain coupling could still bring up unexpected results. Although there 
are two different structural sites of Cu in LiCu202 as shown in figure [T} nonmagnetic 
Zn^+ should substitute the Cu^+ site only from considerations of both valence and ionic 
radius. The ionic radius of Zn^+ ion (0.74 A CN=4) is similar to that of Cu2+ ion (0.71 
A CN=4) for square planar coordination within the Cu02 chain, which is significantly 
larger than that of a Cu+ ion (0.60 A CN=2) within the CuOs dumbbeU.[HI] Room 
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Figure 4. (Color online) Lattice parameters of \j\Cu2-z'^^z02 with various Zn 
substitution levels, (a) The c axis is shown in filled circle and volume is shown in 
filled triangle, and the circle and triangle in (b) represent b- and a-axes, respectively. 



temperature lattice parameters for LiCu2-2Zn^02 as a function of Zn substitution level 
is summarized in figure [4] We note a significant reduction of c axis in concomitant with 
slight increase of a and h axes when crossing the Zn ^ 5 % boundary, with minimal 
change on the cell volume. These results suggest the success of Zn substitution to the 
Cu^+ instead of the smaller Cu+ site within the bridging 0-Cu-O dumbbell, and the 
average lattice distortion must be intimately correlated to the occurrence of novel phase 
transitions found for z > 0.05 at low temperature to be discussed further in the text. 

3.1. Magnetic susceptibilities 

Figure [5] summarizes the complete evolution of Zn substitution effect shown in Xab(T) 
and dXah/dT(T) data for all samples studied. Low Zn substitution reveals typical 
signature of anisotropic helimagnetic ordering transition and reduces slightly below 
Zn ^ 5.5 % per CUO2 chain to show peaks of dx(T)/dT at ^ 20 and 22.5 K, which are 
slightly lower than that of the pristine crystal but similar to the nearly Li deficiency- 
free sample as reported previously. fTO] Across the Zn ^ 5.5 % phase boundary, x(T) 
transforms from a spiral spin ordering character toward a step-like crossover or a cusp 
shape near the phase transition as shown in figure [5](b). The different characters for the 
helical and the novel phase transition are displayed not only by the symmetry change 
of dx(T)/dT peak shape but also by its smooth crossover (< 5.5 %) versus cusp shape 
(> 5.5 %) difference in x(T)- 

The lowering T/^ for helical spin ordering with Zn < 5.5 % per CUO2 chain suggests 
that the helical spin ordering condition of \J2/Ji\ > \ must be met still, which is not 
surprising considering the perturbed spin chains through low level nonmagnetic ion 
substitution. Low Zn substitution does not break the original spiral spin ordering but 
only reduces its transition temperature slightly from ^ 22 K to the lowest ^ 20 K 
before two coexisting phases emerge. It is interesting to note that the single phase 
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Figure 5. (Color online) Temperature dependence of dXah/dT for LiCu2-2Zn202 
along the ah plane with 2; = 0, 1.6, 3.4, 5.8, 8.1, and 10.9%. For clarity, the signal for z = 
10.9% is multiplied by three. x(T)ab data are shown in (b), which are normalized to the 
x(300K) value. The helical ordering temperature versus Zn content is summarized 
in (c). Finite size scaling behavior is found for Zn below ~ 5.5% for the two transition 
temperatures with H applied along the ab plane and c direction. 




Figure 6. (Color online) Magnetic susceptibilities for Zn ~ 5.8 % sample with field 
applied along the ab and c directions. An estimate of susceptibility reduction is 
obtained from Ax/x ^c- 
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sample of the lowest has been achieved through Lio.99Cu202 before, and similar ^ 20 
K transition temperature is found when the Li deficiency is reduced to the minimum. [lOj 
Clearly the spiral spin ordering of LiCu202 system prefers an optimum level of Cu^^ 
excess 0.17/f.u.) instead of exactly equal amount of Cu+ and Cu^+, which must be 
strongly correlated to the required incommensurate periodicity of spiral spin ordering 
with modulation ( = 0.174 along the chain b-direction, while the periodicity is intimately 
related to the J2/J1 ration classically. [9j Current phenomenon of reduction as a 
result of Zn substitution could be due to the finite chain effect, i.e., the Zn substitution 
cuts off the original helically ordered infinite chains before the optimal incommensurate 
modulation for the spiral ordering is broken. While Ji may be frustrated locally near 
the Zn site, the average J2/J1 ratio may not be modified very much at Zn levels lower 
than ^ 5.5 %, which in fact has been supported from our HTSE fitting results to be 
discussed below. 

The Zn substitution above ^ 5.5 % per CUO2 chain introduces an intriguing phase 
transition which is hard to identify based on its x(T) character of nearly isotropic 
cusp shape as shown in figure |6] Since the newly discovered phase transition is nearly 
isotropic for this quasi-two-dimensional system of specific spiral plane, it is unlikely that 
the transition to be antiferromagnetic ordering or even spin glass. The step-like isotropic 
reduction of spin susceptibility implies the existence of a gapped phase that comes from 
part of the system, judging from its non-zero susceptibility at low temperature, even after 
a background subtraction of expected enhancing Curie contribution from Zn introduced 
isolated spins. As shown in figure [6} the ratio of x step reduction versus peak x value 
at Th {Ax/x Tc) is very close to the Zn substitution level between ^ 0.045 - 0.055, 
which implies that only part of the system falls to the gapped ground state and is 
proportional to the Zn level, i.e., high Zn substitution introduces a gapped phase that is 
proximate to the Zn centers. On the other hand, more pronounced hysteresis is observed 
in particular below 5 K, which implies that glassy or irreversible domain effect occurs 
below Tc. Bobroff et al. have proposed that Zn impurities could introduce collective spin 
freezing around each impurity center for spin gapped system at low temperature, [32j 
which may also apply to the gapped phase in our system as suggested by the spin 
glass like behavior below ^ 5 K. Further muon spin resonance experiment is planned to 
examine this proposed explanation. 

We traced magnetic field dependence of transition temperature of crystals with Zn 
level higher than ^ 5.5 % and found the scaling of l-Tc(H)/Tc(0) shows a convincing 
IP dependence as seen in the inset of figure [7| The transition temperatures have been 
determined by linear extrapolation from the upper and lower parts of the susceptibility 
data near Tc before they deviate again due to Curie contribution, similar to that found 
in the Zn substituted spin Peierls compound CuGe03.[33j Bulaevskii et al have shown 
that commensurability effect in spin- Peierls (SP) transition would cause strong magnetic 
field dependence of Tc in the limit of /llbH ^ kBTc(O) to follow equation as[34j 



1 



Tc{H) 



/5[ 



kBTcioy 



(2) 
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Figure 7. (Color online) Magnetic susceptibilities of single crystalline LiCu2-2Zn2 02 
measured near the phase transition temperature Tc(H) in the range of H = 0.01 - 
4.5 T. Tc(H) value is defined by the intersection of two linear fitted lines from above 
and below. A universal plot of l-Tc(H)/Tc(0) versus [-iiBH/kBTc(0)]^ ^ including two 
samples of different Zn levels above ^ 5.5 %, is shown in the inset with a linear fitting 
of slope ~ 0.26. 



Such field dependence is a result of band filling cliange tliat is sensitive to the competing 
commensurability-incommensurability requirement between spin and phonon for spin- 
Peierls transition. We find that there is clear IP dependence for the reduced critical 
temperature. Such If field dependence has been observed in CuGeOa of well known 
spin-Peierls (SP) transition before. [33] The quadratic field dependence found in equation 
[2] also ruled out the possibility of ordinary structural transition due to phonon instability 
under strong field. [34j There is no evidence of structural transformation or signature of 
lattice doubling found based on our synchrotron x-ray structure analysis for the whole 
Zn substituted series (not shown) . The proportional constant /3 shown in equation [2] is 
fitted to be ^ 0.26, which is lower than the 0.46 found in the dimerized system CuGeOa 
but higher than the calculated value of 0.11 within Hartree-Fock approximation. [331 [34] 
Moreover, we find such quadratic field dependence of Tc is independent of Zn content as 
shown in figure [7[ which implies that the same universality class applies to the impact 
of Zn substitution above ^ 5.5 %; and this implication must be clarified further in the 
future both experimentally and theoretically. 



3.2. Finite size scaling 



Figure [5| shows that helical ordering transition temperature Th(z) for LiCu2-2Zn^02 
decreases with higher z before it drops precipitously for z ^ 0.055 critically. The 
precipitous drop of T}^(z) suggests a power law z dependence of large exponent. Fitting 
Th(z) with reduced temperature and z in power law as 
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and we find the exponent n to be near 4 as T^fO)^ Zc and n are left as free parameters 
in the fitting. An equaUy good fitting is obtained when n is fixed at 4, and Th(0) and 
Zc return to be 22.3(3) (for H\\ab), 24.3(3) (for H\\c) and 0.057(2), respectively. The 
satisfactory fitting of equation |3] suggests that finite size eflFect is able to describe the z 
dependence of T^. Finite size scaling has been applied in Monte Carlo simulation using 
limited size L and its size extrapolated exponent u to calculate the critical temperature 
in real infinite size matrix as described in equation |3| [35j Since the critical exponent u 
within mean field approximation to be 1/2, the fitted n=4 suggests that reduction is a 
result of limited correlation length due to finite size, i.e., L(z) oc 1/^. The most natural 
model that can explain such z dependence of L(z) would be a picture of two-dimensional 
helical ordered domains which are cut oS by defect or secondary phases formed domain 
boundaries, the higher the Zn level the smaller the domain size ^ a(b)/2^. Assuming 
dopants are distributed at the two-dimensional domain boundaries, the eflFect of lower 
can be explained as a result of cut-off correlation length by the domain boundaries. Since 
the in-plane magnetic correlation length ^ab for the pure LiCu202 has been estimated 
to be - 1200 A,[9J which is very close to the domain size estimated for z ^ 0.055 at the 
phase boundary, i.e. L = 6/0.055^ ^ 1000 A. It is implied that long range helimagnetic 
spin ordering is destroyed totally once the domain size becomes smaller than the helical 
magnetic correlation length. 

The original undoped sample has a quasi-two-dimensional matrix of helical ordered 
incommensurate spin structure with periodicity about 5.7b, i.e., modulation vector ( ^ 
0.174 along the b-axis with pitch angle 62.6°. [9] We find that the phase boundary for 
Zn ^ 5.5% is very close to the onset of non-zero probability of aggregated (more than 
one) Zn per 6bx2a magnetic supercell of incommensurate helical ordering, i.e., starting 
from 2 out of ^ 35 Cu sites are substituted by Zn. The estimated phase boundary 
agrees very well with the phase boundary Zc obtained from the finite-size scaling fitting 
as discussed previously. The phase boundary obtained from simple probability analysis 
suggests that long-range helical ordering is suppressed only when more than one Zn exist 
per unit superlattice of incommensurate helical ordering. For Zn substituted sample 
prepared from high temperature melt, it is reasonable to assume that Zn ions distribute 
randomly. Although it is unreasonable to assume Zn to be mobile at low temperatures, 
we suspect that the comparable and competing (frustrating) couplings Ji, J2 and J± 
could transport spin degree of freedom from Zn centers to the domain boundaries that 
surround helically ordered spins. Such phenomenon could be viewed as a microscopic 
phase separation with localized spins confined within domain boundaries and encircle 
domains with spins of persisted helical ordering, an interesting contrast to the phase 
separation resulting from one-dimensional charge stripe separated antiferromagnetic 
domains in La2-xSra;Cu04.[36j 
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Figure 8. (Color online) High temperature series expansion (HTSE) fitting and 
exact Hamiltonian diagnolization using 16 spin ring simulation results of magnetic 
susceptibility data for samples with Zn ~ 0, 3.4 and 5.8 %. 



Table 2. Exchange coupling constants Ji and a=J2/Ji obtained from high 
temperature series expansion (HTSE) method and exact diagonalization of the 
Hamiltonian with 16 spin ring (N=16) simulation. Estimations from neutron scattering 
spin wave analysis and LDA calculation are listed for comparison. 



Zn content 


HTSE 


N=16 


Neutron study[14j LDA[IH] 


(EPMA) 


Ji(meV) 


a 


a 


Ji(meV) a Jj_ Ji(meV) a J± 





-6.23(1) 


-0.63(2) 


-0.578 


-5.95^ -0.62t 0.9^ -8 -1.777 5.7 
-7* -0.535* 3.4* 


0.016(3) 


-6.22(2) 


-0.62(1) 






0.034(2) 


-6.21(3) 


-0.61(1) 


-0.54 




0.038(2) 


-6.20(2) 


-0.60(3) 






0.058(2) 


-6.18(2) 


-0.48(4) 


-0.468 




0.081(3) 


-6.12(3) 


-0.42(2) 






0.109(3) 


-6.00(1) 


-0.38(1) 







I model 1 of reference [14j 
model 3 of reference [14j 
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Figure 9. (Color online) Ji, a and J2 (converted from a) obtained from HTSE fitting 
results versus Zn substitution level for 1j\Cm2 -z'^'^z^ 2- 

3.3. HTSE and exact diagonalization fittings 

In order to trace the frustrating couplings Ji and J2 {oi=J2/Ji) for this helimagnetic 
system, we applied high temperature series expansion (HTSE) and exact Hamiltonian 
diagonalization methods to fit the magnetic susceptibility data for the entire Zn 
substitution range. Starting from Hamiltonian written as equation [T} the HTSE uses 
expansion coefficients for the frustrated quantum spin chain model calculated by Biihler 
et a/. [37] Exact diagonalization of the Hamiltonian using 16 spin ring (N=16) calculation 
is used to provide satisfactory temperature dependent fitting at high temperature 
range after peak normalization and background subtraction. Both HTSE and N=16 
fittings return satisfactory results on high temperature range and agree very well on 
the estimation of Ji and a as shown in figure [8] and summarized in figure [9] and 
table [2j These results are in agreement with those obtained from LDA calculation 
and neutron scattering spin wave analysis. p!3l fT4] The fitting results using N=16 finite 
chain are in agreement at high temperature range with HTSE within 6(1) %, although 
the low temperature region deviates significantly as a result of gapped finite size chain 
approximation. We find Ji decreases smoothly, crossing the ^ 5 % Zn phase boundary 
by ^ 4 % only, while J2 (calculated from Ji and a) reduction is far more pronounced. 
The significant J2 reduction through Zn substitution moves this frustrated system closer 
to the quantum critical point near a = -1/4 as shown in figure [3| and complex periodicity 
may have occurred. [231 

While high Zn substitution is proposed to introduce isolated dimers into the original 
long-range helical spin ordered system, the magnetic susceptibility measurement results 
should refiect such assumption, i.e., the final form oi x(T) for Zn > 5.5 % must include 
contributions from a gapped state of isolated dimers as well as the isolated spins that 
shows paramagnetic behavior. We fit the high Zn data additionally with combined 
HTSE, isolated dimer and Curie contributions as shown in figure [T0| following equation 




Figure 10. (Color online) Combined HTSE, isolated dimer and Curie- Weiss model 
fitting. 



X Xo ~^ ^chainXchain ~^ ^ dX Dimer ~^ ^ CwXcW ^ 



= Xo + N^Ha^n}^Y.^n.kC^\A^T + Nn ^^^^^^ + N,^-^.{^) 

The only constraint used in the fitting process is the normalization of 
Nchain+N dimer +Ncurie=^- The complctc fitted rcsults are summarized in table [3| Tem- 
perature dependence of isolated dimers with spin gap size Jd follows the Bleany-Bowers 
equation, which can be reduced to describe the low temperature spin gap behavior 
and the Curie- Weiss behavior for dimers at high temperature. [38j The fitted results 
using equation [4] are not perfect, possibly due to the incomplete HTSE chain descrip- 
tion, neglected interchain coupling, and the possible short-range coupling among dimers; 
however, it does provide a better fitting below for higher Zn above ^ 5.5 %. The 
amount of dimer phase A^^ is nearly proportional to the Zn content level, which supports 
a picture of Zn introduced isolated dimers favorably. Moreover, the fitted values of Jd 
agree very well with the LDA and neutron scattering estimation of (see table [2]), 
which supports a scenario of dimer formation as a result of 0-Cu+-0 bridged nontriv- 
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Table 3. Magnetic susceptibility fitting results of high Zn doping level by the proposed 
model of combined spin chain, isolated dimers and spins following equation. |4] 





/J — W. WUl -L j 


/J — W. WOl -L j 


^ — W. J. J. I J. J 


chain 


0.82 


0.76 


0.52 


Jc (K) 


37.5(2) 


37.0(3) 


29.5(5) 


N 


0.1 


0.14 


0.35 


e(K) 


-31.6(9) 


-25.6(4) 


-49.4(7) 


c 


0.62(6) 


0.62(4) 


0.60(2) 




0.08 


0.10 


0.13 


(K) 


45.8(6) 


48.7(5) 


91.8(3) 


^x/xm (%) 


4.5(1) 


7(1) 


0.5(3) 



ial interchain coupling J^. The Curie contribution from isolated spins, ISIqw/-^ increases 
drastically for z ^ 0.11, which is reflected on the weaker and broadened phase transition 
as shown in figure [5](a) also. Following the same phase separation model suggested here, 
we find that the magnetic susceptibility reduction at (Ax/x(^c)) also reflects the 
spin gapped phase proportion as summarized in table [3| before it is smeared by the 
significantly enhanced Curie contribution at z ^ 0.11. Both the xi"^) model fitting 
and the susceptibility reduction at have demonstrated satisfactory agreement on the 
fraction of isolated dimer phase and the Zn content for z > 5.5 %, which supports the 
interpretation of the special magnetic susceptibility anomalies near (see figure [s]) 
based on a scenario of phase separation. 



34. NMR results 

To gain further insight on the nature of the anomaly observed in the temperature 
dependence of x foi" heavily Zn-doped crystals, we also carried out NMR measurements 
for ''Li (nuclear spin / = 3/2, nuclear gyromagnetic ratio ^7^/27r = 16.546 MHz/Tesla) 
in LiCu2-2Zn^02 {z ^ 0.058) by applying an external magnetic field B along the crystal 
c-axis. We chose the NMR carrier frequency / = 33.095 MHz so that we could compare 
our NMR results for Zn-doped sample with an earlier NMR report for undoped LiCu202 
by Gippius et al [13j At 100 K, we observe a typical NMR lineshape expected for nuclei 
with spin 7 = 3/2 with three distinct peaks split by the c-axis component of the nuclear 



quadrupole interaction, Uq. As shown in figure 11, the central peak near B = 2 Tesla 
arises from the nuclear spin 1^ = +1/2 to —1/2 transition, while two additional peaks 
near B = 1.997 Tesla and 2.003 are from the 1^ = ±3/2 to ±1/2 transitions. From the 
splitting, 0.0030 Tesla, between the peaks, we estimate Uq =^ 7^ • 0.0030 = 0.05 MHz. 

Upon further cooling, all three peaks progressively become broader. In principle, 
there are two possible mechanisms of the NMR line broadening in a paramagnetic 
material, i.e., magnetic effects or quadrupolar effects. In the first scenario, all three 
peaks should broaden by an equal amount, while only the satellite peaks broaden in 
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Figure 11. (Color Online) Representative results of field-swept ^Li NMR lineshapes 
measured at / = 33.095 MHz, (a) at 20 K or above, and (b) at 20 K or below. The 
external magnetic field was applied along the crystal c-axis. Notice that the scales of 
horizontal axis are different by a factor of ten between the two panels. 



the second scenario. The observed lineshapes at 50 K and 30 K clearly show that all 
three peaks are equally broadened, hence we conclude that there exists slight magnetic 
inhomogeneity. This is consistent with the fact that Zn^+ ions substituted Cu^+ sites 
to remove the electron spin 5" = 1/2, thus introducing random disorder in the matrix of 
S = 1/2. The splitting Uq between the peaks remains unchanged within experimental 
uncertainties down at least to 30 K. This implies that EFG (Electric Field Gradient) at 
the ''Li site does not change significantly in this temperature range. Below 20 K, the line 
broadening becomes so severe that we can no longer observe the quadrupole splitting 
between the three peaks. At 10 K and 4.2 K, we observe an extremely broad lineshape 
with a very broad peak centered around 2 Tesla and two equally spaced shoulders near 
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Figure 12. (Color Online) The temperature dependence of (a) the FWHM of the 
^Li NMR lineshape, and (b) the nuclear spin-lattice relaxation rate 1/Ti measured at 
the center of the lineshape. 



1.92 Tesla and 2.08 Tesla. We summarize the temperature dependence of the FWHM 



(FuU Width at Half Maximum) of the NMR hneshape in figure [T2| The FWHM begins 
to increase dramatically below ^ 18 K, where x figure |6] shows a well-defined kink. 

It is very instructive to compare these NMR lineshapes with those observed for 
undoped LiCu202 by Gippius et al fT3] under nearly identical experimental conditions. 
Our results above 30 K are very similar to the undoped case, but our lineshape at 4.2 K 
is strikingly different from that observed for LiCu202. In the latter case, Gippius et 
al demonstrated that two large "horns" develop at 1.95 and 2.05 Tesla, accompanied 
by two additional horns near 1.9 and 2.1 Tesla (see the upper panel of figure 3 by 
Gippius et al). They also demonstrated that these four horns are the consequence of 
the distribution of hyperfine magnetic fields at ''Li sites, arising from an incommensurate 
spiral modulation of magnetic moments. Our NMR lineshape at 4.2 K does not exhibit 
such sharp horn structures, but the overall shapes do bear a strong resemblance to 
each other. The smearing of the sharp horns in our case is naturally understood as 
the consequence of disorder induced by Zn ions in the incommensurate modulation of 
magnetic moments. 

The establishment of this low temperature modulated state with strong disorder 
is unlikely to be due to a second order phase transition. In the case of a conventional 
second order magnetic phase transition, the nuclear spin-lattice relaxation rate 1/Ti 
would diverge at the phase transition temperature because of the critical slowing down 
of the magnetic moments. In contrast, 1/Ti measured at the center of the ''Li NMR 
lineshape does not diverge near 18 K in the present case. Instead, we observe a broad 
hump centered around 12 K. This is consistent with glassy development of the modulated 
ground state. Such glassy behavior has also been demonstrated by the pronounced 
thermal hysteresis in the magnetic susceptibility measurement shown in figure [H] It 
is reasonable to assume that the NMR results describe properties of the major spin 
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Figure 13. (Color online) Specific heat of LiCu2-z'Z^^z02 with Zn = 0, 0.034, and 
0.081. Xab foi" Zn = 0.081 is displayed also to show Tc defined by the cusp shape 
anomaly. 



population in the system, besides the isolated dimers or spins suggested above. The 
glassy development below ^ 12 K could be coming from the frozen short-range ordered 
spins within the frustrated finite spin chains. 

3.5. Specific heat 

Specific heat for samples with different levels of Zn substitution are shown in figure 



13, The undoped and low Zn (< 5.5 %) samples show typical double peaks which 
correspond to the two distinct transition temperatures along the ab- and c-directions 
respectively for the spiral spin ordering. Low Zn sample of 3 % substitution does not 
alter the transition width and simply shift the transition doublet to slightly lower, which 
indicates the robust nature of the helical spin ordering on zero spin perturbation. On 
the other hand, a significantly broader specific heat anomaly is observed near ^ 14 K 
for Zn ~ 8 % sample, in fact there are two very broad peaks sitting above and below 
Tc as defined by the cusp of x(T). These broad anomalies suggest that either the phase 
transition is of short range in nature, or significant inhomogeneity exists in this Zn 
substituted crystal. The existence of broad Cp peaks together with the clear x(T) cusp 
observed at Tc suggests that short-range ordering of spin glass-like behavior occurs. 
Although it is tempting to interpret the coexisting broad Cp anomaly and x(^c) cusp to 
be coming from a well-defined spin glass phase transition, we cannot ignore the fact that 
there are two broad peaks of Cp across Tc, besides, the thermal hysteresis of x(T) is not 
pronounce enough immediately below Tc and only developed at much lower temperature 
as suggested by both the 1/Ti and x(T) below ^ 5-10 K. A reasonable interpretation 
can be obtained following our proposed phase separation scenario, i.e., we might view 
the broad "doublet" to be coming from short-range ordering of the frustrated spins (of 
reduced J2/J1 ratio) on the finite spin chain, which has been cut short by the Zn ions, 
while the cusp at x(^c) is a result of spin gap opening from those isolated spin dimers. 
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Figure 14. (Color online) Electric polarization and dielectric response along the c- 
direction for LiCu2-^Zn^02 with Zn = 0, 0.034, and 0.058. An electric field of 1200 
kV/m has been applied above the transition temperature before each measurement. 



Magnetic entropy ASm can be extrapolated from Cp after the subtraction of 
lattice contribution and integrating from 2 to 160 K. The entropy difference between 
the pure and the ^ 8% Zn substituted samples is estimated to be ^ 580 mJ/mole 
K, which is rather close to a portion (z) of the total Cu^+ spins in the system as 
zR'^ln(2S+l)=0.081^8.3Hn2=470 mJ/mol K, i.e., due to the missing spins for gapped 
singlet formation as a result of Zn substitution. It is puzzling to note that the integrated 
ASm^ 18 J/mole K is almost three times larger than the theoretical expectation, i.e., 
ASm = ^In 2 = 5.8 J /mole K for a 5^=1/2 system. This discrepancy must be due to 
certain factors other than phonon or spin. A possible candidate could be the long- or 
short-range ferroelectric ordering. 



3. 6. Electric polarization 

Dielectric constant and electric polarization have been measured for LiCu2-^Zn^02 {z 
^ 0, 0.034, 0.058) single crystals along the c-direction and shown in figure 14 The 
undoped sample shows the onset of electric polarization below ^ 22 K and two anomalies 
near 22 and 24 K in dielectric constant similar to those reported in the literature. [3j 




Figure 15. (Color online) A proposed dimer configuration for high Zn substituted 
(> 5.5 %) sample below the phase transition temperature. This figure is constructed 
by randomly distributing Zn on a two-dimensional spiral ordered background with 
a incommensurate 2ax6b magnetic superlattice (modulation ( = 0.174) periodicity 
along the b-axis. The alternating dashed orange and solid blue colored lines are not 
in the same plane but bridged by Cu+ as a quasi-two-dimensional bilayer-A as shown 
in figure [l] 



3.4% Zn substitution reduces the spiral ordering temperature and electric polarization 
slightly. 5.8% Zn substitution sample shows no trace amount of electric polarization near 
its magnetic transition temperature near 14 K as indicated in figure [5} Spin-current 
or inverse Dzyaloshinskii-Moriya mechanism has been applied to interpret the origin 
of electric polarization through spiral spin ordering and has been tested by chirality 
reversion experiment. j2[ |3] Current dielectric constant measurement results once again 
confirm that long-range spiral spin ordering is a prerequisite for the spontaneous electric 
polarization in current system. 



3. 7. Proposed dimer models 

For a system composed of coexisting isolated spin dimers on a short range ordering 
helical background, figure [15] shows two simplified dimer models we proposed to 
describe the observed gapped experimental results for Zn substitution level higher than 
5.5%. In the first possible scenario, singlet dimers could form due to non-negligible 
antiferromagnetic interchain coupling (or Jw^ see figure [T]) for interchain coupled 
spins near Zn pairs. EflFective antiferromagnetic interaction between spins near Zn 
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impurity located on different chains can also be nonzero for coupled frustrated quantum 
spin chains, where four-spin coupling (cyclic exchange) has been considered. [391 [40] Since 
the ferromagnetic Ji is disrupted locally by Zn impurity and the antiferromagnetic J2 
becomes dominant locally, the second possible scenario can be described by dimers 
formed with antiferromagnetic coupled Cu^+ ions ( J2) across the Zn impurities within 
the same chain (see figure 15). Such gapped excitation described by a confined spin 1/2 
soliton has various gap sizes as spin chain is cut short from the gapless infinite system to 
various finite lengths. [41j The distribution of gapped excitations can also be refiected on 
the non- vanishing spin gap, broad Cp peaks and ''Li NMR 1/Ti. Isolated spin dimers or 
aggregated clusters coexist on the short-range ordered finite spin chain background when 
Zn level is above ^ 5.5 % below Tc. The proximity of Tc defined hy x(T) reduction onset 
and the broad short range ordering signature in Cp(T)^ as described in figure 13, suggest 
that the coupling strength of short range ordering background spins is comparable or 
even prerequisite to the isolated spin dimerization near Zn pairs. 

For Zn level below ^ 5.5%, dilute isolated spins introduced by Zn should distribute 
randomly near doped Zn centers. However, the observed finite size scaling behavior 
discussed above seems to suggest a scenario that these isolated spins could "transport" to 
the domain boundaries encircled the long-range helical ordering spins, while Zn impurity 
should not be mobile at such low temperatures. Potential spin transport mechanism 
could be described as a soliton, which can be formed near Zn center and transport 
through the frustrating helical ordered spin background. Such behavior supports a 
picture of finite size eflFect for Zn < 5.5% and the short-range ordered dimer cluster 
formation for Zn higher than 5.5%, both of which require spin transport assisted phase 
separation on the frustrating background. Mechanism for creation and propagation 
of a chiral soliton pair on a spiral ordered spin chain background has been described 
theoretically. [H] While nonmagnetic Zn ions break the original infinite spin chain of 
helical spin ordering, soliton pair formed around Zn centers could be transported through 
this intrachain gapless excitation mechanism until it hit the next Zn impurity, which 
could be a starting point on interpreting our findings in this study and a detailed 
calculation is necessary. 

Peierls transition has been one of the most common electronic instability 
mechanisms for the one-dimensional system. There are two types of spin dimer 
formation, with or without spin-phonon interaction, i.e., whether there is phonon 
assisted local structure distortion as singlet dimer forms. Since the phonon assisted spin- 
Peierls transition requires translational symmetry breaking, i.e., the lattice doubling 
which loses translational symmetry to open a gap near the Fermi level, which is a 
collective phenomenon for the one-dimensional spin system. Although the magnetic field 
dependence of shown in figure [7] suggests that the gapped phase has a spin Peierls 
system character, we have not found evidence of lattice doubling from synchrotron 
x-ray diffraction results so far. In fact, it is not expected to be observed based on 
phase separation model, i.e., the phase fraction of isolated dimers is proportional to 
the Zn level of < 10% only. On the other hand, we find that J2 and J± are both 



Nonmagnetic impurity perturbation to the quasi-two-dimensional quantum helimagnet LiCu20223 



antiferromagnetic and of comparable strength according to LDA calculation and neutron 
scattering spin wave analysis. [131 [Hj Spins near Zn pair centers must face the frustrating 
interaction between antiferromagnetic J2 and J±. Interestingly, several LDA and spin 
wave analysis models suggest the ratio of antiferromagnetic J±/J2 to be near 1/2, 
another condition that satisfies the Majumdar-Ghosh dimerization without phonon 
assisted lattice distortion. [42] Whether the dimerization occurs with or without phonon 
interaction, a local probe on subtle structure distortion is necessary; however, it would 
be difficult to verify if only a fraction z) of the spins to demonstrate spin dimerization. 

Microscopic phase separation due to extrinsic inhomogeneity can often been 
observed in samples prepared in powder form or thin film on a substrate. On the 
other hand, phase separation of electronic origin could be an intrinsic property in 
a strongly correlated material system and has attracted studies from many aspects, 
particularly for superconductor and colossal magnetoresistance (CMR) materials. [43j 
Since current finite size effect and spin dimer glass behavior are observed in single crystal 
sample with randomly distributed zero spin ions on a quantum frustrating helimagnetic 
background, the phase separation phenomenon in Zn substituted LiCu202 system must 
be of electronic origin. In particular, even though the driving force remains to be 
explored theoretically, separated phases distinguished by different spin order /disorder 
types are unique and have never been reported before to our knowledge. 

It is interesting to note that the phase diagram summarized in figure [5] shows 
high similarity to phase boundary and Zn dependence for the general phase diagram 
as summarized by Bobroff et al on Zn substituted spin gapped systems. [32j For the 
gapped systems including isolated ladders, Haldane or spin-Peierls chains, Zn impurities 
are proposed to induce extended magnetic moments that freeze at low temperatures. 
Zn impurities must introduce similar residue spins as shown in figure 15 on the spiral 
ordered background, which could explain the spin freezing phenomenon observed near 
^ 5 K from NMR and magnetic susceptibility measurement results discussed above. 

The existence of singlet spin dimer has been predicted from DMRG-type algorithm 
calculation by Furukawa et al for one-dimensional quantum frustrated spin chain with 
a Hamiltonian similar to equation [T| where nontrivial easy plane exchange anisotropy is 
considered also in order to simulate the asymmetric ''Li NMR spectra accurately. [131 US] 
We note that the high Zn substituted samples with J2/J1 ratios summarized in table [2] 
could fall in the range of dimerization according to the phase diagram constructed by 
easy-plane anisotropy vs. J1/J2 as proposed by Furukawa et al. In particular, the J2/J1 
ratio drops significantly as system crosses over the ^5.5% Zn boundary from the original 
chiral spin ordering to the partially gapped spin dimerization state. However, we must 
point out that the phase diagram proposed by Furukawa et al is based on an infinite 
chain system calculation while current Zn substituted spin chains are finite of unequal 
lengths, which could explain why the dimerization does not open the spin gap fully but 
with a size proportional to the Zn substitution level. 

Finally, the current results also bring in an interesting connection to the possibility 
of experimental realization of the long sought resonant valence bond (RVB) state. For a 
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zigzag spin chain system, it has been shown that the RVB state is the exact ground 
state at the ratio of J2/Ji= -1/4. [47j Since figure [9] shows that Ji remains nearly 
invariant while J2 decreases with the increasing Zn, one might view our system is 
actually approaching this RVB critical point at high Zn level, i.e., Zn level higher than 
^5.5% is near the onset of the RVB state. The RVB state on a finite spin chain due 
to spinless Zn substitution has gapped low-lying excitations which give rise to a cusp 
shaped susceptibility as indicated in figure [5)3. Clearly more theoretical exploration is 
necessary when interchain coupling and easy-plane anisotropy are considered. 

4. Conclusions 

In conclusion, we have used nonmagnetic Zn substitution to explore the quantum quasi- 
two-dimensional helimagnet LiCu202. Zn substitution moves the frustrated spin system 
closer to the quantum critical ferromagnetic boundary near a = -1/4: on the J2 vs. Ji 
phase diagram. Two kinds of magnetic orderings have been found crossing the T-z 
phase boundary for z near ^ 5.5 % per CUO2 chain. Low Zn substitution reduces 
the spiral ordering transition temperatures significantly, following a finite size scaling 
law, which implies the existence of two-dimensional helimagnetic domains confined by 
isolated spins, and these isolated spins could be transported as a soliton through the 
frustrating spin background to the domain boundaries. Such spin phase separation 
is uniquely observed for the first time in a two-dimensional helimagnetic background. 
High Zn substitution introduces a gapped phase that shows spin-Peierls character and 
that the phase fraction is proportional to the Zn substitution level. The proposed 
interpretations on phase separation, spin dimerization and spin freezing require more 
microscopic experimental evidences and muon spin resonance investigation of the title 
compound is under way. 
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